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SUMMARY 

The growth and decomposition rates of nuclei on a surface are expressed 
In microscopic terms. A stability factor Is defined and employed for the 
derivation of a criterion for nucleus stability. Simulation results Indicate 
that the stability factor can be used as a measure of the system stability. 


INTRODUCTION 

The nucleatlon of thin films has been the subject of Intensive research 
for a long time (ref. 1). For some reason, however, only a small fraction of 
this effort was devoted to the problem of nucleus stability, a central topic 
In classical nucleatlon theories (refs. 2 and 3). Moreover, some of the few 
studies which do address this problem (refs. 4 and 5) are based on the classi- 
cal theories (refs. 2 and 3) and utilize macroscopic quantities such as the 
free energy and surface tension. As pointed out by Walton (ref. 6), however, 
the applicability of those theories to nucleatlon of thin films Is questiona- 
ble, since experimental observations (refs. 7 and 8) have Indicated that the 
critical nucleus In this case Is small and consists of but a few atoms. It Is 
thus evident that a different approach, microscopic in nature, Is needed In 
order to treat this problem. Such an approach Is presented here. 


THEORY 

The change In a nucleus size Is the outcome of two competing processes, 
namely growth caused by joining of atoms, and decomposition due to detachment 
of atoms. We deal first with the growth process. The kinetics of this proc- 
ess depends on the rate of arrival of atoms to the nucleus boundary. When an 
atom joins the nucleus, a bond (or several bonds) Is formed. Thus, this proc- 
ess results In a decrease In energy. Following Metropolis et al . (ref. 9) we 
assign a value of one to the joining probability. In other words, we assume 
that every atom which arrives at the nucleus boundary joins the nucleus. When 
the flux J of atoms arriving at the nucleus boundary Is homogeneous, the 
rate A of arrival of atoms to the nucleus Is given In macroscopic terms by 

A = JL (1) 

where L Is the boundary length. 

In microscopic terms the equivalent of L Is the number of available growth 
sites, l.e., the number of accessible lattice sites Immediately adjacent to 
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the nucleus which we designate by ng.. Thus, the growth rate 6 Is given 
In microscopic terms by 


G = 3n g (2) 

We now examine the decomposition process. The jump frequency of an 
adsorbed atom on the substrate surface Is given by 

» ■ v xp (jf) < 3 > 

where t> 0 Is a characteristic frequency, E the binding energy, k Is 
Boltzmann's constant, and T the absolute temperature. In the nearest- 
neighbor model (ref. 10), which we employ here, the binding energy E Is 
assumed to be the sum of single bond energies, l.e.. 


E 



+ Z 


(4) 


where E a Is the energy of binding between the atom and -the substrate, E& 
the energy of a single bond between adsorbed atoms, and Z the number of 
nearest neighbors. Thus, 


v 


v 0 exp 



+ Z 

Tr 



(5) 


or 



where given by 


( 6 ) 


v d 


v 0 exp 
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Is the frequency of jumps of single atoms (l.e., those which are not bound In 
clusters) . 

The flux of atoms detaching form the nucleus, l.e., the decomposition rate D, 
Is therefore given by 


0 . (8 

where the sum Is over all the nucleus atoms located on the boundary. 

The tendency of the nucleus to Increase In size can be expressed by means 
of the ratio between the growth and decomposition rates: 


G 

D 



(9) 


2 



When this ratio Is greater than unity the nucleus has a tendency to Increase 
In size and vice versa. 


The quantities J and vjj In equation (9) depend on the conditions on 
the substrate and on the nature of the growing material, respectively. The 
quantity which depends on the nucleus Is the ratio which we will henceforth 
call the stability factor and designate by s: 



Equation (9) can now be written In the following concise form 



( 10 ) 
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The flux J Is the produce of the concentration and the velocity. Expressing 
It In terms of the number of atoms crossing a unit of length equal to the lat- 
tice parameter In a time unit one obtains 

J = 0V(j (12) 

where e Is the surface concentration. 

Substitution to equation (11) yields 



(13) 


A nucleus Is considered stable when Its tendency to grow exceeds Its 
tendency to decrease In size. We thus obtain the following criterion for 
nucleus stability 


es > 1 (14) 

We now demonstrate the use of this criterion In two simple examples of 
nuclei on a two-dimensional square lattice. 


Linear Nuclei 

For a linear nucleus composed of n monomers and lying along the (1,0) 
direction the number of available growth sites Is 2n + 2, and the stability 
factor Is 


2n + 2 

s = 2 

(n - 2) cl + 2a 

where a Is the Boltzmann factor 



(15) 
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Thus, the stability criterion In this case Is 

(n - 2)a 2 ± 2a 
6 2n + 2 

This form of the criterion Is valid, of course, only for the Initial configu- 
ration. Once an atom joins the nucleus or detaches from It, the nucleus shape 
changes and a different form of the criterion has to be used. 


Square Nuclei 

For a square nucleus composed of n^ monomers and bound by four (1,0) 
lines the number of available growth sites Is 4n, and the stability factor Is 


2 3 

a + (n - 2)a 

Thus, the stability criterion In this case Is 

e > a 2 + (n - 2)a 3 


(18) 


(19) 


Finally, It should be emphasized that throughout this treatment repre- 
sents the local surface concentration In the Immediate vicinity of the nucleus. 
This Is equal to the nominal coverage only before extensive growth has taken 
place. At later stages depletion zones are formed around the growing nuclei, 
and this Is no longer true. 


COMPUTER SIMULATION 


The above equations can be easily utilized for computer simulation. By 
choosing l/u^ In equation (6) as the time unitone obtains the following 
expression for the probability P of each atom (excluding those which are 
completely surrounded by nearest neighbors) jumping at any time 



( 20 ) 


Simulation studies using this jump probability were conducted on a square 
lattice containing 10,000 (100 by 100) sites. Eight simulation runs were per- 
formed for each data point and the average as well as the standard deviation 
were calculated. 


Since we are Interested now In the stability of the system as a whole, we 
define a system stability factor S In a manner similar to that used to define 
the nucleus stability factor, namely 



( 21 ) 
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where f Is the coordination number of the lattice (In our case f = 4) and 
the sum Is over the whole system. 

The variation of this factor with time for a surface coverage of 0.1 and 
several values of Ejj/kT Is presented In figure 1. A feature which Is common 
to all of these curves Is that for all values of Eb/kT the stability factor 
Is a nondecreasing function of time. This Is In contrast to other parameters 
such as the percentage of bound atoms or the cluster density which were found 
(ref. 11) to decrease with time when Eb/kT < 0.5. 

Comparison of these curves with the resulting microstructures, presented 
In figure 2, shows that once S levels off (e.g., for Eb/kT = 2 at t = 100) 
no detectable change In microstructure takes place. It may thus be concluded 
that the factor S can be used as a measure of the system stability. 


CONCLUSIONS 

The stability of nuclei on a substrate surface has been expressed In 
microscopic terms, namely the binding energy and the local atomic configura- 
tion. A stability factor has been defined and a criterion for nucleus stabil- 
ity has been derived. The results of computer simulation Indicate that the 
stability factor can be used as a measure of the system stability. 
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Figure 1.- Variation of the system stability factor S 
WITH TIME FOR SEVERAL VALUES OF E r /KT. 
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Figure 2.- Evolution of microstructure with time for several 
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